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In the context of a simple microscopic schematic scalar model we study the effects of spatial 
correlations in force transmission in granular assemblies. We show that the parameters of the 
normalized weights distribution function, P(v) ~ v a exp(— v/(j>), strongly depend on the spatial 
extensions, £y, of such correlations. In particular we evaluate the functions <j>(£v) an d a(£v)- If 
£v — » oo then <j) ~ £v and weights are power law distributed. 



Photoelastic visualizations of stress propagation in granular media have provided evidence of strong inhomogeneities 
■ in the distribution of contact forces with the formation of "force chains" which typically extend on space scales much 
larger than the dimensions of single grains These inhomogeneities results from local fluctuations and are 

accompanied by broad stress distributions observed in bead packs which are responsible for many unusual 

I | properties of granular media |jj . 
i Several theoretical approaches have been proposed to describe the unusual properties of stress patterns in such 
systems [[|-[7j and recently a schematic scalar model, which takes into account just the vertical component of the 
t-H , stress tensor, has been introduced to describe the statistics of stress distribution pL||. The mean field study of this 
model has shown rich results in good correspondence with experimental observations. However recent experiments ||] 
show the necessity to go beyond the mean field approximation and to take into account spatial correlations present 
in force transmission from grain to grain, correlations which actually lead to chains formation. 

In the present paper, we study the effects of such correlations in the framework of a schematic scalar model, as 
in ref. [§,[§, in which grains, disposed on a lattice, discharge their weights in random fractions, due to the disorder 
' present in a typical granular assembly, on their bottom neighbours. The disorder of the environment, however, does 
Q\ , not imply that force transmission between grains is uncorrelated |^,p|,^2|. In the complex dynamical process which 
' leads an assembly to a static configuration, a certain fraction of sites carrying a strong shear (typically involved in 
the formation of "force chains"), slips away from some of their own neighbours, losing contact with them. These 
processes may induce strong spatial correlations in force transmission, which we try to take into account here. 

In our model the weight, w(i, h), supported by a grain at height h and column i of a square lattice (see Fig. [I]) is 
(using the notation of ref. [|]||): 



w(i, h) = q + (i — 1, h — 1) w(i — 1, h — 1) 

+q-(i + l,h-l)w(i + l,h-l) + l (1) 



where q+ (resp. g_) is the fraction of the weight that site (i — 1, h— 1) (resp. (i + l,h — 1)) discharges on site (i,h), and 
we have supposed that the masses of single grains are equal to unity. Due to the conservation of the mass, we have 
the constraint q+(i, h) + q~(i, h) — 1. The values of q+(i, h) are generally uniformly distributed in the interval [0, 1]. 
However, we suppose that a fraction 1 — 8 (5 G [0, 1]) of the total number of grains (a fraction of grains randomly 
displaced on the lattice) is subjected to the slip condition proposed in ref. ||: 

q+(i,h) = 0(x) (2) 

Here 9 is the Heaviside step function and x is |Io| : 

x = [q+(i — 1, h — 1) w(i — 1, h — 1) 

-q-(i + 1, h - 1) w(i + l,h- 1)] /w(i, h) (3) 

Eq. (|^) schematically expresses the slip condition in which, due to friction, a grain, strongly "pressed" from its top 
right (resp. left) neighbour, may discharge its weight mainly on its own right (resp. left) lower neighbour pl] |. This 
mechanism originates force chains in our model. 

Our tuning parameter is the fraction of grains not subject to the slip condition, 8, a parameter which, in real 
samples, is related to friction coefficients, elastic constants, grains shapes. The model interpolates the two extreme 
situations corresponding to a complete uncorrelated weight transmission, studied in ref. Q, which is obtained when 
5 = 1, and the strongly correlated case studied in ref. ||, which is recovered in the limit 5 = 0. 
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This two dimensional model is simple enough to visualize mechanisms underlying force transmission. Actually, the 
sites which undergo the above "slip condition" are typical points enhancing correlations in force transmission and so 
belonging to "force" chains (see Fig. |]) . In our case the "slipping" grains are randomly distributed on the lattice; thus 
the probability to have a segment of such grains of total length d, i.e., to have a chain of length d, is approximately 
V{d) ~ (1 — 5) d . Consequently the average chain length, £y, which is experimentally easier to measure, is related to 
the fraction of such sites as: £y ~ (1 — 8)/5. The kind of correlation we have introduced above essentially generates 
"chains" which sustain the weight of the lattice columns which intersect them from above, as "arches" HJj. Thus, in 
the deep bulk of the system at depth L, the weight, W, that a point at the base of a chain of length d experiences, is 
approximately W ~ Ld. A quantity of theoretical as well as experimental and practical importance is the distribution 
P(v) of normalized weights, v = W/L, felt by grains at a given depth L. With the above arguments we can derive an 
approximate expression to relate P(v) to the fraction 5, i.e., to the average chain length £y: 

P(v) = P(W(d)/L) = V(d) ~ exp(-v/<t>) (4) 

where <f> ~ l/log[l/(l — 6)] {4> ~ £y, if £y — ► oo). This result implies important practical consequences, because it 
predicts that we have to expect to measure weight fluctuations of the order of £y . Thus if £y , which is a characteristic 
length of the specific system we consider, is very high, huge stress fluctuations arise in the sample. The above picture 
compares well with the numerical calculations we present below, where its consequences are discussed in more details. 

Our numerical analysis of the present model concerns a square lattice of length L = 10000 and depth M = 1000, 
averaged over 10 different realizations. In this lattice we adopted periodic boundary conditions along the horizontal 
axis. 

To study the effects of spatial correlations introduced by the slip condition of eq. (^) , we evaluate the vertical space 
correlation function of forces, Cy(r), along a main axis of our square lattice: 

(w(i, M)w(i -r,M - r)) - w m (0)w m (r) 
Lv[r) " {w^Mf)-w m {Q) 2 [ ) 

with w m (r) = {w(i,M — r)}. The behaviour of Cy(r), in the present model, is depicted in Fig. | as a function of 
the distance between grains r (expressed in lattice units). After a first jump from Cy(0) = 1 to Cy(l) ~ 0.5, Cy(r) 
smoothly decreases approaching zero. The first part of the decay is exponential in r, and may be fitted with the 
following function: 

C v (r) =K v cxp(-r/Z v ) (6) 

The characteristic length £y of eq. (|^) , is related to the measure of the extensions of forces inhomogeneities along the 
vertical direction in our system. It is the typical vertical extension of forces correlation, or, more crudely, the average 
vertical distance between crossing points of "stress paths" observed in photoelastic measurements (see Fig. [TJ). In 
Fig. ||, £v is plotted as a function of 5 [[u}. In agreement with the above theoretical arguments, the approximate 
behaviour £y ~ <5~ 1,0 is found, showing that £y increases orders of magnitude when the fraction of sites, 1—6, 
undergoing the slip condition approaches 1. In the limit £y — » oo, force chains extends over the whole system and 
this fact strongly affects weights distribution. 

The horizontal space correlation function, Cjj(r), is defined as the average over the lower 10% of the system (to 
have more precise data) of c(r, h) : 

(w(i,h)w(i + r,h)) - (w(i,h)) 2 
(w{i,hy) - (w{i,h)y 

The function Ch{t), depicted in Fig. ^, becomes negative as soon as r > 1, signaling that horizontally aligned grains 
are always slightly anti-correlated. When r — > oo, Ch(t) asymptotically approaches zero from below. As above, the 
first part of the decay is almost exponential: 

C H (r) = -K H exp(-r/S H ) (8) 

In this case, the length corresponds to the typical horizontal spacing of chains. The function £h{£v) 1S reported 
in Fig. |3|. diverges when S — > 0, and we approximately find £h ~ £y with a ~ 1/2. 

We now try to relate the above observations to the distribution, P(v), of the weights, v = w/w m , normalized by 
the mean w m — h, supported by grains at a given depth h. The function P(v), evaluated at the bottom layers in our 
model, is plotted in Fig. ^. In agreement with the above theoretical considerations, this quantity becomes independent 
of h if measured at sufficient depth in the system. In ref. has been proposed that P(v) is characterised, in a very 
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broad class of models, by the following behaviour, which is approximately recovered in the present model (we find 
deviations at very small v): 

P(v) = Av a exp(-v/^)) (9) 

An interesting result from our numerical calculation, consistent with the theoretical considerations presented above, 
is that the asymptotic exponential behaviour of P(v) is recovered for all values of 5 > 0. The correlations present in 
force transmission between grains, except exceptional cases, do not alter this property. The validity of this important 
observation may be broader, being consistent with results obtained in different contexts as the carbon paper experiment 
with compressed stationary bead packs of ref. [0] or experiments on continuously sheared granular materials of ref. [0 , 
or in numerical simulations of scalar or vectorial force models [| 12 l^JIT 



However, the quantities a and <j), which in mean field theory seem to be exclusively related to the coordination 
number of the lattice of grain packing |§|,|l2]| , strongly depend on the degree of correlation present in forces transmission, 
i.e., in our model, on the fraction S. This fact has important practical consequences and may explain contrasting 
results found in different experiments. The dependence of the distribution P(v) on spatial correlations has been 
recently experimentally outlined in ref. S. 

In the present model, the exponent a decreases in presence of spatial correlations in force transmission between 
grains. It passes from the value a = 1 at S = 1 (as predicted by mean field theory) to the value a ~ — 1.1 at S = 
(in agreement with the simulations of ref. ||). At a fraction of "normal" sites S ~ 0.7, a value higher than the 
percolation threshold on the square lattice, a crosses the zero. The exponent a is depicted in Fig. |as a function of 
the vertical correlation length £y . Figure |5| shows that in the region of not too large values of £y , small changes in £y 
may induce large variations of a. Such a strong dependence seems to be found also in the experimental observations 
of ref. H , where was shown the sensitivity of the measure of a with respect to the relative sizes of the grains and 
of the measuring device. In particular, from the above general discussion and the results of ref. we expect that 
the results of mean field theory have to be experimentally recovered essentially when the measures of forces are taken 
averaging over regions which are larger than the length £y. 

For what concerns the parameter <fi, Fig. [s| shows that it diverges approximately as a power law of £y The 
superimposed fit, consistent with the previous theoretical arguments, is </>(£v) ~ Thus 4> grows with £y and in 
the limit of huge spatial correlation, i.e., £y — > oo (or S — > 0), the exponential asymptotic decay of eq. (|^) is lost and 
just the power law behaviour survives. This fact has significant practical importance because it implies that giant 
stress fluctuations can be observed 0. 

In conclusion, we have studied a simple microscopic model in the framework of scalar force approximation 
which schematically takes into account the effects of correlations present in forces transmission in granular assemblies. 
Its simplicity allows to clarify the effects induced by the extension of "forces chains" in the system. The presence of such 
a non trivial characteristic length scale, which is experimentally measured for instance by photoelastic visualization, 
can thus be quantitatively related to other measurable quantities as the forces distributions, P(v). We have discussed 
some interesting correspondences with known results from recent experiments . However further experimental and 
theoretical investigation about these important effects is still missing. 

Interestingly the present model is related (see ref. ||) to works devoted to describe other physical phenomena as 
directed Abelian sandpiles, aggregation-dissociation reactions, interface dynamics, or river networks Eq]. 
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FIG. 1. Two pictures of the conformation of stress paths in our model, corresponding to two values of the correlation length, 
£v (£v ~ 2 in the left picture, and £y ~ 20 on the right). They are obtained enlightening the grains which carry a weight 
above a given threshold. 

FIG. 2. The vertical, CV(r), and the horizontal, Cfr(r), force correlation functions as a function of the grain distance 
r (in unit of the lattice spacing), for several values of the fraction of grains not undergoing the "slip condition" of eq. (Q), 
8 — 0.1, 0.3, 0.6, 1.0 (resp. squares, triangles, diamonds, circles). The superimposed curves are the exponential fits described in 
the texts. 

FIG. 3. Left: the characteristic length £y of the vertical force correlation function as a function of 8. It diverges when 
<5^0 approximately as ~ 5~ 10 . In such a limit force correlations extend over the whole system, and this fact strongly 
affects forces distribution. Right: the characteristic length £h of the horizontal force correlation function as a function of £y 

(€h~#"). 



FIG. 4. The distribution P(v) of forces v = w/w m normalized by the mean force w m at the bottom of our system, depicted 
for three values of the fraction of grains not undergoing the "slip condition" of eq. (^), 8 = 0.1, 0.6, 1.0 (resp. squares, diamonds, 
circles). The superimposed fits are from eq. (g). 

FIG. 5. The parameters of eq. (^), a and (f>, as a function of the correlation length £v- Right: the exponent a passes 
from the value predicted by mean field theory, a = 1, at small £y (i.e., 8 = 1) to a — —1.1 when £v — > oo (i.e., at 8 — 0). 
The sensitivity of a to changes of £v remembers the observations from experiments by Miller et al. [3]. Left: the parameter 
<f> diverges as a power law with £y (approx. <f) ~ showing that if £y — > oo the exponential asymptotic decay of force 
distribution P(v) is lost, and huge stress fluctuations are possible. 
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